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We discuss the limits of validity of the semiclassical theory of gravity in which a 
classical metric is coupled to the expectation value of the stress tensor. It is argued 
that this theory is a good approximation only when the fluctuations in the stress 
tensor are small. We calculate a dimensionless measure of these fluctuations for a 
scalar field on a flat background in particular cases, including squeezed states and 
the Casimir vacuum state. It is found that the fluctuations are small for states which 
are close to a coherent state, which describes classical behavior, but tend to be large 
otherwise. We find in all cases studied that the energy density fluctuations are large 
whenever the local energy density is negative. This is taken to mean that the gravita- 
tional field of a system with negative energy density, such as the Casimir vacuum, is 
not described by a fixed classical metric but is undergoing large metric fluctuations. 
We propose an operational scheme by which one can describe a fluctuating gravita- 
tional field in terms of the statistical behavior of test particles. For this purpose we 
obtain an equation of the form of the Langevin equation used to describe Brownian 
motion. 
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I. INTRODUCTION 



A natural proposal to describe the gravitational field of a quantum system is the semi- 
classical theory in which the expectation value of the stress tensor is the source. The 
semiclassical Einstein equation is 

Gfiu = Rfiu — -^diivR = ^T^GNiTfMu) ■ (1-1) 

This theory is almost certain to fail at the Planck scale, where the quantum nature of 
gravity becomes important. However, it can also fail far away from the Planck scale if the 
fluctuations in the stress tensor become important. This was discussed some time ago by 
one of us and illustrated by the problem of graviton emission by a box of particles in 
a general quantum state. It was found that the semiclassical theory gives reliable results 
when the fluctuations in the stress tensor are not too large, that is, when 

{T^p{x)T^M) ^ {T^p{x)){T^M)- (1-2) 

However, for quantum states in which the energy density fluctuations are large, the semi- 
classical theory based upon Eq.( |LlD cannot be trusted. 

In the present paper, we will further explore the issue of the limits of validity of the 
semiclassical theory. Particular attention will be paid to quantum states in which the ex- 
pectation value of the local energy density can be negative. In Sec. y some aspects of such 
states will be reviewed, partly for the purpose of establishing notation and results which will 
be used later. In Sec. |T| the range of validity of the semiclassical theory will be probed in 
particular flat space examples, including squeezed states and the Casimir effect. In Sec. |^ 
a proposal will be made as to how to describe the gravitational field of a system in which 
the energy density fluctuations are large. We suggest that a statistical description of the 
Brownian-like motion of test particles in such a gravitational field yields all of the informa- 
tion which is available about the system. Finally, in Sec. the results will be summarized 
and discussed. 

We adopt the convention c = h = 1 and a spacelike metric rj^y = diag[—l, 1, 1, 1] 
throughout this paper. 

II. QUANTUM STATES OF NEGATIVE ENERGY DENSITY 

Although all known forms of classical matter have non-negative energy density, it is not 
so in quantum field theory. A general state can be a superposition of number eigenstates and 
may have a negative expectation value of energy density in certain spacetime regions due to 
coherence effects, thus violating the weak energy condition 0. If there were no constraints 
on the extent of the violation of the weak energy condition, several dramatic and disturbing 
effects could arise. These include the breakdown of the second law of thermodynamics 0, 
of cosmic censorship and of causality citeMorris88. There are, however, two possible 
reasons as to why these effects will not actually arise. The first is the existence of constraints 
on the magnitude and the spatial or temporal extent of the negative energy [0,^0. The 
second is that the semiclassical theory of gravity may not be applicable to systems in which 
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the energy density is negative. This latter possibihty will be the main topic to be investigated 
in this paper. 

Let us consider a massless, minimally coupled scalar field for which the Lagrangian 
density is 

C^\r^,.{d<^ct>){d^ct>). (2.1) 

The stress tensor is 

T^, = n(5o0) - 7]^.C = {d^(t>){dA) - \riAda4>W(t>), (2.2) 

where IT = d^cj) is the conjugate momentum of 0. Variation with respect to (f) gives the 
dynamical equation 

(-^ + V2)0(a:) = d^d^<t>{x) = 0. (2.3) 
The quantum field operator cf) may be expanded in mode functions as 

= EK/k + aU), (2.4) 

k 

where 

[aic, Ok/] = 5kk', [ok, Ok'] = [al, al,] = 0, (2.5) 

and the mode function is a solution of the dynamical equation 

/k = {2L^u)-^/''e"'^''^ = (2L^a;)-^/'e''^-"-''^*. (2.6) 

Here we have assumed periodic boundary conditions in a 3-dimensional box of side L. The 
sum J2k is transformed to / (Pk in the continuum limit. The mode functions are nor- 
malized on a spacelike hypersurface under the scalar product 

(01, 02) = -^ / r-'x{M^)docl>;{x) - (ao0i(x))0;(x)}, (2.7) 

so that 

(/k,/k') = <^kk'. (2.8) 

We will here consider states where a single mode is excited. Such a quantum state takes 
the form 

00 

l*) = Ecn|n), (2.9) 

n=0 

where \n) is a number eigenstate with n particles in mode k and the c„ arc any coefficients 
such that Xl^o = 1- The normal ordered expectation value of the stress tensor is 
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oo 

= ^ (2n|cjX/3[/k, /k] + n'/'in - lY^'cnC:_,To^p[h, /k] 

n=0 

- l)^/^c:c„_2T.;,[/^, /^]), (2.10) 
where ^^/^[si, /i] is the bihnear form 

r^4^,/i](x) = - ^r^^.(9,^)(9'^/i). (2.11) 

The different bihnear forms arising from the mode functions are 

Tapih^M^-ICo^pe''', (2.12) 

^a/3[/k) /k] = Ta/3[fk, /k] = •'Cq,/?, (2-13) 

and 

T^pif^j:] = -fC^fse-'"' , (2.14) 

where ^ = kpX^ and /Cq/3 = (/Cq-Zc/^ — ^rjapkpk'') / {2ujL^) . For the massless case fepfc'' = 0, 

so }Caf3 = (^xf)' Noti'^^ that :Tai3:= T^p — (0|Ta/3|0); the normal ordered stress tensor 
in flat spacetime is the renormalized result obtained by subtracting the Minkowski vacuum 
expectation value. 

The quantum coherence effects which produce negative energy densities can be easily 
illustrated by the state composed of two particle number eigenstates 

I*) ^ ^7^^(10) + e|2)), (2.15) 

where |0) is the vacuum state satisfying a|0) = 0, and |2) = ;^(a''")^|0) is the two particle 
state. Here we take e, the relative amplitude of the two states, to be real for simplicity. For 
this state, 

{■.T^p{x):)= {%:T^p{x):\m) 



e 



{ V2(r,^[/k, /k] + i;M/k, /k]) + 2e(r,^[/k, /k] + i;M/k, m 



l + e2 

= -^^(2e - v^cos(2^)). (2.16) 
1 + e"' 

Obviously the energy density can be positive or negative depending on the value of e and 
the spacetime-dependent phase 9 = kpX''. We also observe that the negative contribution 
comes from the cross term. For a general state which is a linear combination of particle 
number eigenstates, the number of cross terms will increase as A^^. Therefore, for a general 
quantum state, the occurrence of a negative energy density is very probable. However, the 
total energy is always non-negative. We will come back to this simple illustrative example 
when discussing the breakdown of the semiclassical gravity theory later. 
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III. LIMITS OF THE SEMICLASSICAL THEORY 



To obtain a criterion of the validity of semiclassical gravity theory, we may recall the 
investigation in Ref. [Q. There the energy flux of gravitational radiation in linearized gravity 
produced by a matter field was calculated both in the semiclassical theory and in a linear 
quantum gravity theory. In the semiclassical theory based upon Eq.( |l.l| ), the flux depends 
on products of expectation values of stress tensor operators, whereas in a theory in which 
the metric perturbations are quantized, it depends upon the corresponding products of 
expectation values. For the particular case of a single-mode coherent state, the product 
of expectation values is the same as the expectation value of products of stress tensors. 
This can be understood from the fact that a coherent state corresponds to a classical field 
excitation, for which one would expect the semiclassical theory to be a good approximation. 

We propose that the extent to which the semiclassical approximation is violated can be 
measured by the dimensionless quantity 



This quantity is a dimensionless measure of the stress tensor fluctuations. (Note that it is 
not a tensor, but rather the ratio of tensor components.) If its components are always small 
compared to unity, then these fluctuations are small and we expect the semiclassical theory 
to hold. However, the numerous components and the dependence upon two spacetime points 
make this a rather cumbersome object to study. For simplicity, we will concentrate upon 
the coincidence limit, x — > y, of the purely temporal component of the above quantity, that 
is 



The local energy density fluctuations are small when A <^ 1, which we take to be a measure 
of the validity of the semiclassical theory. Note that we have used normal ordering with 
respect to the Minkowski vacuum state to define the various operators. 

It is not difficult to see why the semiclassical gravity theory is not expected to be valid 
when the energy fluctuations are large. Suppose we have a quantum state which is a super- 
position of two states, each of which describes a distinct macroscopic matter configuration, 
e.g., (1) the presence of a 1000 kg mass, or (2) the absence of this mass. Clearly, a mea- 
surement of the gravitational field should reveal either (1) the field of a 1000 kg mass, or 
(2) no field, each with 50% probability. But the semiclassical field equations predict finding 
the field of a 500 kg mass with 100% probability. 

Now we will discuss some specific examples involving the massless scalar field in fiat 
spacetime. First we consider the general quantum state l^') = Z^^o Cn\n), for which a single 
mode is excited. The correlation function for the stress tensor can be divided into three 
parts. 



(:T„^(x)r^^(y):) 



(3.1) 




(3.2) 



( 



(3.3) 



The contribution of the diagonal terms is 
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Po = Y.\<^n?{n-lf''n{n + lf/' 

n=0 

(^a/3[/k) /k] Tfj,u[f^, /k] + Taplf^-, /k] ^Mi'l/k) /k] 
+^a/3[/k, /k] ^M^^i/k) /k] + ^a/3[/k, /k] ^A'l'l/k, /k] 

+T„^[/k, /k] T^.[/k, /k] + T«M/k, /k] ^,..[/k, /k]). (3.4) 
That of the matrix elements between \n) and |n + 2) is 

oo 

^2 = E 4+2 n (n + 1)^/2 (n + 2)^2 

n=0 

(T'a/ji/k, /k] ^Ati't/k' /k] + ^Q/3[/k5 /k] T^y[f^, 

+T^^[f:, T^4/k, /k] + r«/3[/k, /k] ^^/..[/k, /k]) + c.c. , (3.5) 
whereas that between |n) and |n + 4) is 

oo 

P4 = E Cn c„+4 {n + 4)V2(^ + 3)V2 + 2)^/^ (n + 1) t„^[/*, /*] t^^[/*, /*] + c.c. , (3.6) 

ri=0 

where c.c. denotes the complex conjugate. The effect of the normal ordering has been to 
remove the contribution of the n = term from the diagonal part, Pq. 



A. The Vacuum + Two Particle State 



Now we go back to the simple case of a state which is a superposition of the vacuum and 
a two particle eigenstate, Eq. (|2.15|) . In this case we find 

(: T„^(x) T^,(x): ) = ^^^^ (T„^[/*, /*] T,,[h, /k] + /J r^.[/k, /k] 

+Tal3[fk, /k] Tf^u[fk, /k] + Tal3[h, /k] T^u[fki /k] 
+^Q/3[/k, /k] ^M'^l/k, /k] + ^a/3[/k, /k] ^/xi^f/k) /k]) 
126^ 

From this result, Eq. ( p.l6| ), and the fact that /Coo = 2l3' ^® ^^"^ ^^^^ 

A(.) = + (3.8) 

From Eq.( p.l6[) , it follows that the condition for the expectation value of energy density to 
be negative is cos {29) > y/2e. In this case we have A{x) > 1. This indicates when the 
energy density is negative, the energy density fluctuations are large and the semiclassical 
theory is not a good approximation. 
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B. Squeezed States 



Squeezed states of light have been extensively investigated recently in quantum optics 
and have been experimentally realized 0. Squeezing reduces the quantum uncertainty in 
one variable with a corresponding increase in that of its conjugate variable. Squeezed states 
form a two parameter family of quantum states, which include both coherent states and 
states with negative energy density as different limits. 

A general squeezed state for a single mode can be expressed as 

\a,O = D{a)S{C)\0), (3.9) 
where D{a) is the displacement operator 

D{a) = exp{aa^ - a*a) = e'l^l'/^ ^aat ^-a*a ^g^^g^ 
and S{Q is the squeeze operator 

SiO^expi^Ca'-lcia^]. (3.11) 

Here 

a = se'^ (3.12) 

and 

C = re'^. (3.13) 
are arbitrary complex numbers. The displacement and squeeze operators satisfy the relations 

D^{a)aD{a) = a + a, (3.14) 

D\a)a^ D{a) = + a\ (3.15) 

(C) a 5(C) = a cosh r - ah'^ sinh r, (3.16) 

and 

S\Oa^ S{0 = a^ cosh r-ae-^"^ sinh r. (3.17) 

Using the above formulae, the expectation value of the renormalized stress tensor is easily 
obtained: 

{a, C\:T^p{x): \a, C) = (0| S\() D\a) (aX/3[/k, /k] + a) a {T^^[f^, /*] 

+T^p[f:, /k]) + (a^)X/.[/k, m) D{a) S{() |0) 
= 2/CQ,/3{sinh r cosh r cos(26' + 5) + sinh^ r 

+s^[l -cos 2(^ + 7)]} (3.18) 
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Simlarly, the expectation value of the squared stress tensor is 

(a, C|: T„^(x)r^,(x): |a, C)= /C^, (^/[cos4(^ + 7) - 4 cos2(^ + 7) + 3] 

+35^ {2 sinh r cosh r [2 cos(2^ + 5 + 27) - cos(4e + 5 + 27)] 
+4 sinh^ r (cos 27 — cos 26) — cos(5 + 27) | 

+3sinhV[coshVcos(40 + 2(5) + 3-4cos2^]y (3.19) 

The degree of squeezing may be measured by the squeeze parameter \(\ = r. When 
C = 0, the squeezed states reduce to coherent states. In this particular case, the expectation 
value of a product of stress tensors is equal to the corresponding product of expectation 
values 

(: T^p{x) T^M- ) = {■■ Tap{x): ) (: T^Av)- ), (3-20) 

and hence 

A„^^,(x,y) = A = 0. (3.21) 

This result is consistent with the interpretation of coherent states as the quantum states 
which describe classical field excitations. By the criterion which we have adopted, the 
semiclassical gravity theory is a good approximation for coherent states. 

The opposite limit from a coherent state is the case where a = 0, known as a squeezed 
vacuum state. Such a state is not, of course, the vacuum state so long as C 7^ 0, but rather 
a superposition of states containing even numbers of particles. Squeezed vacuum states 
have a particular physical interest because they are the states which result from quantum 
particle creation processes. Such states of the electromagnetic field have been generated in 
the laboratory using nonlinear optical media, and have been the topic of much interest in 
quantum optics in recent years. 

In a squeezed vacuum, a = 0, we may also take 5 = 0, as this is simply a choice of phase, 
and write 

{a,(\-Tai3{x): \a,() = 2/Ca/3 sinh r [cosh r cos(2^) + sinhr]. (3.22) 

Here r = corresponds to the vacuum state and hence gives a vanishing expectation value for 
the stress tensor. For r 7^ 0, the squeezed vacuum state exhibits negative energy densities. 
That is, for fixed r as ^ varies from to 2tt (either the spatial position changing by one 
wavelength at fixed time, or the time varying through one period at fixed position), the 
energy density becomes negative during part of the cycle. 

For the case 6 = and a = 0, the expectation value of the squared stress tensor is 

(a, C|: Tai3{x)T^iy{x): \a,() = 2/Ca/3 ICfj,^ sinh^ r[2 cosh^ r cos46' — 8 sinh r cosh r cos 26 

+3(sinhV + coshV)]. (3.23) 

As required, for the vacuum state (r = 0) this quantity vanishes. 
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From the above results, we can form the quantity A. However, the analytical expressions 
are not particularly transparent. Our primary concern is whether or not A ^ 1, which is 
best determined by numerical evaluation of A using Eqs. (|3.2|), ( |3.18| ), and ( |3.19| ). The 
figures illustrate the results. For states which are sufficiently close to coherent states, i.e. 
r -C I a I, we do indeed find that A ^ 1. However, as the magnitude of the squeeze parameter 
r increases relative to that of we find that A increases. By the point that the state is 
sufficiently squeezed to have p = {: Too :) < 0, we always have that A is at least of order 
unity. Thus squeezed states for which the energy density is negative exhibit large energy 
density fluctuations. 



C. The Casimir Vacuum 

One of the most astonishing predictions of quantum field theory is the Casimir effect 
0, in which the vacuum energy density of the quantized electromagnetic field between two 
parallel perfectly conducting plates is negative. We will consider the Casimir effect for a 
massless, minimally coupled scalar field. The expectation value of the stress tensor can be 
expressed as 

{T^u{x)) = I lim (V^. V. + V^V.. - ^7^.V,V"') G{x, x'), (3.24) 

where the Hadamard elementary function G{x, x') is 

G{x,x') = ^{(f){x)(p{x') +(p{x')(p{x)). (3.25) 

This quantity is formally infinite and is renormalized by replacing G{x, x') by the renormal- 
ized Green's function 

Gr{x, x') = G{x, x') - Go{x, x'), (3.26) 

where Go{x, x') is the Minkowski space Green's function, i.e., that in the absence of bound- 
aries. The resulting finite stress energy is the expectation value in the Casimir vacuum of 
Tfj^i, normal ordered with respect to the Minkowski vacuum: 

(: T^,{x) :)c = ^ lim (V^. V, + V^V,. - ^7^,V,V"') ^^(x, a;')- (3.27) 

Z x'^x 

For the purpose of illustration, we will consider the particular boundary condition of 
periodicity in the 2;-direction with periodicity length L. In this case, the Green's function 
can be expressed as an image sum: 

oo 

G{x,x')= GnMx.x')), (3.28) 

n=— oo 

where 

G„,(a„(x,x')) = -^. (3.29) 
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The geodesic distance a for different n is 

a„(a;,a;') = -(^-^')' + (r-r'J^ (3.30) 

where 

X = (t, r), x' = it', r'), = (x', y\ z - riL). (3.31) 

The renormahzed Green's function Gr{ OC ^ X ^ IS simply given by Eq. (|3.28|) with the n = 
term omitted. 
Thus we find 

(: T^u{x): )c= ^ lim i^t^'^^^ + ^i^^'^' + ^t^u^^''9pda')a~'^ 

n=— oo 

C(4) 

- ^ diag[l, -1,-1,3], (3.32) 



where (^(s) = S-^i is the Riemann zeta function, and we have excluded the n = term, 
since it is just the contribution from the Minkowski vacuum. Since ({A) = 

(:Too(a:):) = -^. (3.33) 

We now proceed to define the stress tensor correlation function as before by interpreting 
■.Taf3{x)T^^{y): as being normal ordered with respect to the Minkowski vacuum state. Its 
expectation value is obtained by differentiation of the corresponding four-point function 

(: T^p{x)T^^{y):) = Jim^^ [VaxV^y - ^Vaf^Vxi) iVt^pVua - ^Vi^uVp^t) x 

di dl dl (: 0(a;i)0(x2)0(x3)0(x4) :)c. (3.34) 

This four-point function is the expectation value of an operator normal ordered with respect 
to the Minkowski vacuum, but whose expectation value is taken in the Casimir vacuum 
state. We have in fact dealt with similar quantities in our calculation of the mean stress 
tensor, namely (:T^;/: )c itself and the renormahzed Green's function 

Gr{x,x') = 0(x)0(x') + 0(x')0(x) :)c. (3.35) 

The evaluation of the expectation value of the four-point function can be achieved with 
the aid of Wick's theorem. Let 0j = 0(xi). The quartic product may be expressed as 

01020304 = : 01020304 : + : 0102 : (0304)m+ : 0103 : (0204)m+ : 0104 : (0203)m + 
: 0203 : (0104) A/+ : 0204 : (0i03)a/+ : 0304 : (0i02)m + 

(0102)A/(0304)Af + (0103) Af (0204) Af + (0104) A/ (0203) Af , (3.36) 

where () u is the expectation value in the Minkowski vacuum. However, it may equally well 
be written as 
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01020304 = iVc(0l020304) + iVc(0102) (0304)c + -^c(0103) (0204)c + 

-^C (0104) (0203) C + ^C(0203)(0104)c + 
-/Vc (0204) (0103) C + ^C(0304)(0102)c + 

(0102)c(0304)c + (0103)c(0204)c + (0104)c(0203)c, (3.37) 

where A^c denotes normal ordering with respect to the Casimir vacuum. Recall that Wick's 
theorem applies to any quantum state |^) for which there is a decomposition of the field 

operator into positive and negative frequency parts, = 0'^+0~, so that = (V''|0~ = 0. 

Both the Minkowski and Casimir vacua satisfy this condition. We now wish to equate the 
two expressions above, solve for : 01020304 :, and take its expectation value in the Casimir 
vacuum. If we use such relations as (-/Vc(0i02))c = and 

(0102)m = (0102)C - (A^m(0102))c, (3.38) 

we obtain the result 

(: 0(a;i)0(x2)0(x3)0(a;4) :)c = (: 0(a;i)0(a;2): )c {■ (pix3)(l){xi): )c + 
(: 0(xi)0(x3): )c {■■ 0(2^2)0(^4): )c + (: (t>{x^)(i>{x^): )c {■ 0(^2)0(^3): )c- (3.39) 

Using this result, we may express the expectation value of the squared energy density as 

(: Tqq^{x): )c = hm {5qx (^07 - T.Vxi) (<^op ^oa - 7:Vpa) d] x 

((: 0(a;i)0(a;2): ) (: 0(a;3)0(a;4): ) + (: 0(a;i)0(a;3): ) (: 0(a;2)0(a;4): ) + 
(: 0(a;i)0(a;4): ) (: (I){x2)(j){x3):)). 



= C7^oo:)^+^[(:0^:)^ + (:(0,.)^:)?; + 

(:(0,,)^:)^ + (:(0,,)^:)^]. (3.40) 

Let p = (: Too: )c be the Casimir energy density and let Pi = (: Tif )c be the pressure in the 
i-direction. Then for the massless scalar field we have 

(:0'Oc = ^(P + Pi+P2+P3)' (3.41) 



and, for example. 



{■■{(p,.f--?c = \{p + Pi-P2-'P,f. (3.42) 



Prom these relations, we may write 

(P + 6 - 6 - 6)' + (p - 6 + 6 - 6)' + (p - Ci - 6 + 6)% (3.43) 

where = pi/ p- The quantity A' is related to A by 
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A A' 

A' = r> A = (3.44) 

1-A' 1 + A'' ^ ' 

and is an equivalent measure of the scale of the energy density fluctuations. 

We may obtain lower bounds on A' and A. The minimum value of the right hand side 
of Eq. (|3.43|) is 1/2 and occurs when = ^2 = ^3 = 0. Thus we have that 

A' > ^, A > i. (3.45) 

This means that the dimensionless Casimir energy density fluctuations are always at least 
of order unity. Note that this bound is independent of the details of the geometry of the 
boundaries. For the particular case of periodicity in one spatial direction, we have from Eq. 
( ^.32| ) that ^1 = ^2 = —1 and ^3 = 3, so that A' = 6 and A = 6/7. The essential point 
here is that neither measure of the fluctuations is small, and hence we conclude that the 
gravitational field due to Casimir energy is not described by a fixed classical metric. That 
the Casimir force is a fluctuating force, and that Casimir's calculation yields the mean value 
of that force, have been emphasized by Barton pUj] . 



IV. METRIC FLUCTUATIONS AND TEST PARTICLES 

The properties of spacetime can be probed by test particles which follow geodesies in a 
classical gravitational field. When the gravitational field is described by a classical metric, 
we can use test particles as an operational probe of the geometry. Our problem is now to give 
a meaning to the gravitational field of a fiuctuating source, such as the Casimir vacuum. We 
propose to replace the description in which the test particles follow fixed trajectories by a 
statistical description in which one only attempts to compute average quantities such as the 
mean squared velocity of an ensemble of test particles. This is the approach which is used 
to describe Brownian motion by means of a Langevin equation. The basic idea is that the 
particle is subjected both to a classical force and to a fiuctuating force. For nonrelativistic 
motion, its equation of motion may be written as 

m^^ = F,(x)+F(x), (4.1) 

where m is the test particle mass, Fc is the classical force, and F is the fiuctuating force. 
The solution of this equation is 



v(t) = v(to) + - fl^cit') + ¥{t')] dt'. (4.2) 



We assume that the fiuctuating force averages to zero, (F) = 0, but that quantities quadratic 
in F do not. Thus the mean squared velocity, averaged over an ensemble of test particles is, 
for the case that v(to) = 0, 

(v2) = ^ /* dt, f dt2 [F,(ti) Fe(t2) + (F(ti) Fit,))]. (4.3) 

Jto Jto 



12 



Typically, the correlation function for the fluctuating force vanishes for times separated by 
much more than some correlation time, tc- In this case, the contribution of the fluctuating 
force to (v^) grows linearly in time. 

We wish to consider the motion of a test particle in a weak gravitational field, so 

Qui^ = Vfiu + V' (4-4) 

where <^ 1. We further assume that the particle's motion is non-relativistic. However, 
we make no assumptions concerning the relative magnitudes of the stress tensor components, 
in contrast to the usual Newtonian limit where one assumes that Too is large compared to 
all other components. The proper time interval along the particle's world line is 

dr^— g^i, dx'^ dx" — goo dt^ + Qij dx^ dx-' 

= dt^ {goo + 2goi v' + gij v' v^), (4.5) 

where v"^ = dx'^/dt is the velocity of the test particle. The dynamics of the test particle is 
determined by the variational principle 

,dT 



/(-^)* = 0. (4.6) 



Therefore, the Lagrangian of the test particle is 

dr 
dt 

= (-1 + hoo + 2hoiv' + v^ + hij v' v'y/\ (4.7) 
Up to terms linear in /i^,^, the Lagrangian is 

L^-l + ^hoo + hoi v' + ^v^ + ^hij v' v^. (4.8) 
The curvature tensor is related to hf^u through 

where h — h°'a- In the Lorentz-Hilbert gauge, 

V'' = (V - Ivi^u h^a)'" = 0, (4.10) 
hf^i, can be determined by the gravitational source through the gravitational field equations 

dpd' V = -167rGiv7;. = -IGttGn (T^u - ^^M-T'^a)- (4.11) 
If we separate the spatial and temporal components, we have 

dpd' hoo = -SttGn {Too + T\) , (4. 12) 
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Opd^hoi^-UnGNToi, (4.13) 

and 

dpd^ hij = -IdnGj, {Tij + ^SijToo - ^^i.-T^). (4.14) 

Define generalized potentials, $, and as linear combinations of the metric pertur- 

bations so that 

/ioo = 2$ + 2*, (4.15) 

hoi = Ci, (4.16) 



The potentials satisfy 



Kj^^ij-2^5ij + 2-^5ij. (4.17) 

dpd" $ = -AttGnToo, (4.18a) 

a^a'' * = -AnGNT'i, (4.18b) 

9^9''Ci = -167rG'Arroi, (4.18c) 



and 



dpd^Cij^-i^^GNTij. (4.18d) 
The Lagrangian of the test particle is 

L^-l + {^ + ^) + Civ' + ^v^ (1 - 2$ + 2*) + (4-19) 

The Euler-Lagrange equation derived from this Lagrangian, after discarding terms higher 
than linear in the velocity or its first derivatives, is 

+ di^ + diiCi v^) - a + I'i (1 - 2$ + 2*) + (-2$ + 2^) + 4- + Cij v' = 0. (4.20) 

We can now solve for the time derivative of velocity and notice that in the weak field 
limit 

(1 - 2$ + 2^ + Cii)"' 1 + 2$ - 2^ - ^ij, (4.21) 

so 

Vi = Ci- di^ - di^ + (2$ % - 2^ Sij - diCj - 4) Vj - Cj{di Vj). (4.22) 
To lowest order 

14 



V^ = Ci- J (d^^x)) dt- J (di^ix)) dt. 



(4.23) 



Equations ( |4.18a[ - [4.18d| ) may be solved in terms of the retarded Green's function 

1 



to yield 



GrpAx) 



47r|x| 



5(|x-t|), 



$(t,x) = Gn d^x 



,3 fToo{t- |x-x'|,x') 



X — X' 



(4.24) 



(4.25a) 



^(t,x) = Gn d-^x' 



,3 , T%(t-|x-x'|,xO 



X — X' 



(4.25b) 



0(t,x) =4G'jv / d^'x' 



Toi{t- 




X — x' 


,x') 




X 


-x' 





and 



3^/ Tj jit |x X I, X ) 

|x — x'l 



^ij(t,x) = AGn J d-x 
The expression for the velocity, Eq. (|]2|), now becomes 
Vi{t)= AGn I d^x' 



Toi{t- 




X — x' 


,x') 




X 


-x' 





-Gat a / / dt 



Tooit' - 


X — x' 


,x')+n(t'- 


X — x' 


,x') 




X — x' 





The mean square of the i-component of the velocity is then 

{V,\t)) = lQCo^O^{t,t)+^^ f dt'dt"C0000{t\t") 

f dt'dt"C,\\t', t") - 8di f dt'Comit, t') 
-Sdi y* dt'Coi/{t, t') - 8d^ y* dt'dt"Cooj^{t', t"), 

where 

C^.af3{t',t")=Gl f d^x'd'x^ 



(4.25c) 



(4.25d) 



(4.26) 



(:T,,(t'- 


X — 


x' 


,x')T,^(t"- 


|x - x" 


,x"):) 




X 


-x'l 


|x-x"| 





(4.27) 



(4.28) 



Note that in Eq. ( |4.27| ), the index i is not summed, but all other repeated indices are 
summed. 

This expression is rather cumbersome to evaluate explicitly, but we can make order of 
magnitude estimates. Note that it is of the general form of Eq. ( |4.3| ), where the correlation 
function of the fluctuating force is of order 
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{FF) ^ m'^C, 



(4.29) 



where C is a typical component of Cijki- If is the correlation time, then we obtain the 
estimate 

{Vi\t)) ^ Ct,t. (4.30) 

To be specific, let us consider a test particle moving in the fluctuating gravitational field 
produced by the Casimir vacuum. We assume that the particle is initially shot parallel 
to and midway between two plates. Classically, the particle will continue to move midway 
between the plates, as the gravitational force vanishes by symmetry. However, it is subjected 
to a fluctuating force which tends to cause it to drift toward one side or the other. Since 
the energy density is of the order of the characteristic force on the test particle is of 

order GNfn/L^, so C ~ Gf^/L^. Furthermore, the correlation time is of order L, the only 
length scale in the problem. Therefore, 




(4.31) 



where Lp is the Planck length. Thus the transverse velocity of the test particle undergoes 
a random walk, with the rms transverse velocity growing as the square root of time. We 
have of course restricted our attention to a particular initial condition for the test particles, 
but a similar analysis could be performed for any other initial condition. The resulting 
probabilistic statements about such quantities as (fj^(t)) sum up all that we can know 
about a fluctuating gravitational field, such as that whose source is the Casimir vacuum 
energy. 



V. DISCUSSION AND CONCLUSIONS 

In this paper, we have argued that the semiclassical theory of gravitation based upon 



Eq.( |l.l| ) is an approximation which assumes that the fluctuations in T^^, are small. This 
assumption can fail far away from the Planck scale. The approximation is valid for sys- 
tems described by coherent states, but breaks down for more general quantum states. In 
particular, states of a quantized field in which the mean energy density is negative seem to 
exhibit large energy density fluctuations. This was illustrated in our model calculations for 
squeezed states and for the Casimir vacuum state. Thus the gravitational field of such a 
system is not described by a fixed classical metric, but rather by a fluctuating metric. 

Because the operational significance of a gravitational field lies in its effects upon test 
particles, the meaning of a fluctuating gravitational field is given by a statistical description 
of ensembles of test particles. This notion was developed in Sec. where we proposed a 
Langevin-type equation to determine such quantities as the mean squared velocity of the 
test particles. 

The analysis in this paper was restricted to a flat spacetime background. In this case, the 
expectation values of quartic operators which appear in the dimensionless measure of the 



fluctuations. A, or in the Langevin equation. Eg. ( 4.27 ), may be defined in terms of normal 



ordered products. The generalization of this analysis to curved spacetime backgrounds will 
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require a renormalization procedure for quartic operator products in such spacetimes, which 
has not yet been developed. Note that an alternative approach to fluctuations is to discuss 
only averaged quantities. If we were to average flelds over flnite space or time intervals, then 
we could deflne a quantity similar to A without invoking normal ordering. This approach is 
used by Barton |T0[ in his treatment of the fluctuations of the Casimir force. In the present 
context, this procedure suffers from the defect that it introduces an arbitrary length or time 
scale into the problem. If, for example, one wished to obtain an equation of the form of 
Eq.( [4.3| ) using averaged rather than normal-ordered quantities, there would be an inherent 
ambiguity in the resulting equation. 

A topic of considerable current interest is the possibility of causality violation by a 
Lorentzian wormhole. Morris, Thorne, and Yurtsever have shown how such a wormhole 
might in principle be constructed using the Casimir vacuum energy. A violation of the 
weak energy condition, and hence negative energy densities, is essential for the existence 
of a wormhole. Hochberg and Kephart W^ have suggested squeezed states as a possible 



source for the negative energy, but an explicit wormhole model using squeezed states has 
not yet been constructed. All of the literature on wormhole models assumes the semiclassical 
gravity theory. However, in light of the considerations in the present paper, it seems to be 
important to include the possible effects of metric fluctuations. Whether they are sufficiently 
large to prevent the creation of a wormhole or its use to violate causality is not clear. Another 
question of interest is whether a breakdown of the semiclassical theory provides an alternative 
explanation as to why negative energy fluxes cannot lead to violations of cosmic censorship 
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FIGURES 



FIG. 1. The energy density p is plotted as a function of the squeeze parameter r and of s = |a|. 
Here 7 = 5 = and 9 = tt/2. When s >> r, the state is close to a coherent state, and /? > 0. 
However, when s << r, the energy density is negative. 

FIG. 2. Here A is plotted for the same range and choices of parameters as in Fig. |l]. Note 
that A << 1 only when s » r, the region of classical behavior. Otherwise, A ~ 1. 

FIG. 3. The energy density p is plotted as a function of the squeeze parameter r and of 7, 
which is the phase of a. Here 6 = 0, s = ^, and 9 = 7t/2. 

FIG. 4. Here A is plotted for the same range and choices of parameters as in Fig. ^. Again 
A << 1 only when r << 1. 
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